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1 Introduction. 



In a previous paper we took the first step towards constructing a renormahzable effective 
potential for the local composite operator (LCO) ^4^ = A'J^A'^'^ in linear covariant gauges. It 
was shown within the algebraic renormalization formalism [21 01 that Aj^ is multiplicatively 
renormalizable to all orders in perturbation theory. At the same time, the anomalous dimension 
of yl^ was explicitly computed to 2-loops in the MS scheme as a function of the gauge fixing 
parameter, a, where a = corresponds to the Landau gauge. The computation exploited the 
fact that in linear covariant gauges the operator does not mix, for example, with ghost operators 
of dimension two with the same quantum numbers. 

The operator has recently received widespread interest in Yang-Mills theory in the Landau 
gauge. Its relevance has been advocated both from a theoretical point of view as well as from 
lattice simulations 11100111111 • Analytic results in favour of a non-zero value for the condensate 
<^^^^ in the Landau gauge have been obtained recently, [Hl^l. Further, the inclusion of quarks 
has been considered in jllj . Motivated by the result of ^^1 it has been shown in ^^1 that Aj^ is 
multiplicatively renormalizable to all orders in the Landau gauge, but its anomalous dimension 
is given by a combination of the gauge beta function, /3(a), and the anomalous dimension, 7^(0); 
of the gluon field, according to the relation jl2j 

lA<a) = - + 7A(a)) , « = ll^ • 

An important consequence of the formation of the {^^^ condensate in the Landau gauge is 
the dynamical generation of a gluon mass mgiuon ~ 500MeV [HI. Lattice simulations of the 
SU{2) gluon propagator in the Landau gauge report m^i^^^ w 600MeV, ^11. Gluon masses 
have also been extracted from lattice methods in the Laplacian [151 [IB] and Maximal Abelian 
gauges jl7j I18j . Earlier the pairing of gluons was discussed in connection with mass generation 
as a result of the instability of the perturbative Yang-Mills vacuum jl9l I2U1 OTt .22. .23. . A 
dynamical gluon mass might also be important, for example, in connection with the glueball 
spectra, [23[ini- A dimension two gluon condensate (Af) was already introduced in [^51, where 
the Coulomb gauge was considered. Furthermore, a dynamical gluon mass is part of a certain 
criterion for confinement introduced by Kugo and Ojima, ^26. .27 . For a recent review see |28j . 

It is no coincidence that the Landau gauge is employed in the search for a gluon condensate of 
mass dimension two. As is well known, there does not exist a local, gauge-invariant operator 
of mass dimension two in Yang-Mills theories. However, a non-local gauge invariant dimen- 
sion two operator can be constructed by minimizing Aj^ along each gauge orbit j29j I3nj I31j . 

^min = {VT)~^ mmu J d'^x ^A^^ where VT is the space time volume and C/ is a generic 
SU{N) transformation. This operator A^-^^ is related to the so-called fundamental modular 
region (FMR), which is the set of absolute minima of {VT)~^ J d'^x(^A^^ . In the Landau 
gauge, d^A"-'^ = 0, so that A^^^ and Aj^ coincide within the FMR. As such, a gauge invariant 
meaning can indeed be attached to (^Aj^^ in the Landau gauge, as it was also expressed in [HI. 

Another interesting property of the Landau gauge is that the operator Aj^ is BRST invariant 
on-shell. If one considered alternative gauges to the Landau gauge, one could search for a class 
of gauges in which the operator ^4^ can be generalized to a mass dimension two operator while 
maintaining the on-shell BRST invariance. Doing so, one should consider a class of non-linear 
covariant gauges, which are the so-called Curci-Ferrari gauges [S^llSSl, where Aj^ is generalized 

to the mixed gluon-ghost operator (^A'^A'^"' + ac°'c°'] [S^ESl. The latter operator is indeed 
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BRST invariant on-shell |,S4[ l,S5j . and has been proven to be multiplicatively renormalizable 
to all orders jSHI, and to give rise to a dynamical gluon mass in the Curci- Ferrari gauge j37j . 
Moreover, in the physical meaning of (^^A'^A^"' + ac°c°^ was discussed, based on on-shell 
BRST invariance. 

In the Maximal Abelian gauge, which is a renormalizable gauge in the continuum |39[ I4(jj . 
one should consider the operator (^^A^A'^^ +^c^c^^ where the group index /3 labels the ofF- 
diagonal generators of SU{N) with /? = !,... ,N{N — 1) and ^ is the gauge parameter of the 
Maximal Abelian gauge. This operator also enjoys the property of being both BRST invariant 
on-shell [341 135j and multiplicatively renormalizable to all orders |36l I41j . Although the effective 
potential for the condensate (^^A^Ai^^ + ^c^^c^'^ has not yet been obtained, we expect it to have 
a non-vanishing vacuum expectation value, which would result in a dynamical mass for the off- 
diagonal gluons. This is based on the close similarity between the Maximal Abelian gauge and 
the Curci- Ferrari gauge and hence the results of '37| . 

More commonly, the Landau gauge is a special case of the well known linear covariant gauges. 
Although the operator Aj^ is not even BRST invariant on-shell in these gauges, it is still renor- 
malizable to any order in perturbation theory This is due to the fact that, thanks to the 
additional shift symmetry, c — > c -|- const, of the antighost in the linear covariant gauges, the 
composite operator A"^ does not mix into the dimension two ghost operator cc. In this article, 

we will construct the effective potential for the dimension two condensate ^^^^ in linear gauges 

and show that a non-vanishing value of (^^fij is energetically favourable, resulting in dynamical 
gluon mass generation. 

The paper is organized as follows. In section 2, we briefly review the local composite operators 
formalism and explicitly calculate the 1-loop effective potential. In section 3, we discuss the 
gauge parameter independence of the vacuum energy which requires an extension of the LCO 
formalism. The behaviour of the gluon propagator is discussed briefly in section 4 whilst we 
provide concluding comments in section 5. 



2 LCO formalism and effective potential for Aj^. 

2.1 Construction of a renormalizable effective action for Af. 

We begin with the Yang-Mills action in linear covariant gauges 

S = Sym + Sgf+fp (2) 

where 

j^ab^Q^^ab_gfabc^c^^ (3) 

is the covariant derivate in the adjoint representation. In order to study the local composite 
operator we introduce it into the action by means of a BRST doublet |2J of external sources 
(J, A), namely 

Sj = sj d\ + |aj) = j d'x Q j4 + XAld^c- + I , (4) 

where s denotes the BRST nilpotent operator acting as 

sA^ = -Dfc\ 
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sc" = 


1 h h 


= 


b\ 


sb" = 


0, 


sX = 


J, 


sJ = 


. 



(5) 

According to the local composite operator technique [HI 1421 I44j . the dimensionless parameter 
C is needed to account for the divergences present in the vacuum Green function (^A'j^{x)Al{y)'j , 
which turn out to be proportional to J^. As is apparent from the expressions © and Q, the 
action {Sym + Sgf+fp + Sj) is BRST invariant 

s {Sym + Sgf+fp + Sj) = . (6) 

As was shown in T, the action {Sym + Sgf+fp + Sj) enjoys the property of being multiplica- 
tively renormalizable to all orders of perturbation theory. 

To obtain the effective potential, we set the source A to zero and consider the renormalized 
generating functional 

exp-iW{J) = J[Dip]eicpiS{J) , (7) 

with 

S{J) = Sym + Sgf+fp + Sct + j d^x I Z2 + (C + ^C) y J ' (8) 

where (p denotes the relevant fields and Sgt is the usual counterterm contribution. Also, 5C, 
is the counterterm accounting for the divergences proportional to J^. The bare quantities are 
given by P 

9o = ZgQ , Uo = ZAOi , Co = ^cC ' '^o = ZjJ , 

(9) 

where Z^C = C, + 5C, and Zj = The functional W(J) obeys the renormalization group 

equation (RGE) 

^"11 + ^^^'^^ + "^"^^'^^"^^'^^'^ / d'^Jjj+^i9\C)^^mJ)=^, (10) 



where 



r,{g\0 = /.|-C. (11) 
From the bare Lagrangian, we infer that 

CoJl = ^^-'{C + K)J\ (12) 



d 
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where we will use dimensional regularization throughout with the convention that d = 4 — e. 
Hence 

m|;C = Vig', = 27A<g'K + 5{g\a) , (13) 

with 

%^a) = (^e + 2-fA'^ig^,a)-P{g^)-^-aja{9^a)—^ • (14) 

Now, we are faced with the problem of the hitherto arbitrary parameter As explained in 
1321 ESI 115 5 setting C = would give rise to an inhomogeneous RGE for W(J) 

M|; + /?(5')^ + «7a(5')^-7A2(5') j d^xJ-^^W{J) = 5{g\a) J d^x^ , (15) 

and a non- linear RGE for the effective action F for the composite operator This problem 
can be overcome by making ( a function of g^ and a so that, if g^ runs according to (3{g'^) and a 
according to ^a{g'^), Cid'^i'^) will run according to ((T^ . This is accomplished by setting C equal 
to the solution of the differential equation 

d ' ^A. 

da. 



/3(5')7r^ +a7a(5',«)— ) C(5',«) = 27^2 (^2)^(52, a) + a) . (16) 



Since ^(3^,0) now automatically runs according to its RGE, W(J) obeys the homogeneous 
renormalization group equation 

The final step in the formal construction of the effective potential for ^^^^ is the removal of the 
j2 terms from the Lagrangian by means of a renormalized Hubbard-Stratonovich transformation. 
By this procedure, the energy interpretation of the effective action is made explicit again and 
the conventional IPI machinery applies. We insert unity written as 

l = ^/[0<^lexp[./A(-y^(^-z4-Z,ayj] , (18) 

with A'^ the appropriate normalization factor, in ^ to arrive at the Lagrangian 

/:(A,,a) = -^F,l + + + > (19) 

while 

exp-iW(J) = J[D^]expiSa{J) , (20) 

SaiJ) = j d^x{c{A^,a) + J^ . (21) 

Now, the source J appears as a linear source term for ^. From and (|20)) . one has the following 
identification 



5W{J) 



' 42 \ / 
J=o \ ^ \5/ 



6J 

where we will not write the renormalization factors from now on. This equation states that the 
gauge condensate ^^^^ is related to the expectation value of the field a, evaluated with the new 
action, / d'^xC{Af^, a), of ((TU]) . 

Although we have not considered the contribution from (massless) quark fields in the previous 
analysis, it can be checked that the results remain unchanged if matter fields are included. 
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2.2 Explicit calculation of the 1-loop effective potential. 



Firstly, we will determine the renormalization group function d{g'^,a) as defined in dJ. All 
the following results will be within the MS scheme. The value for P{g'^) can be found in the 
literature. In d dimensions, one has 

Pig^) = -£5' - 2 (/3o5" + + 0(5' 

Pi = (f C^l - 4CpT^iV/ - f C^T^iV/) . (23) 

where the Casimirs of the colour group are defined by Tr(T"'T^) = Tpd""^, j^arpa _ (j^j^ 
jacdjbcd _ (j^^ab g^^^j^ jg ^^iq number of quark flavours. For 70(5^), we use the relation 
la{g^) = ~'^'lA{g^)- The anomalous dimension 7^1(5^) of the gluon field in linear covariant 
gauges was calculated at three loops in MS in [l^l ■ Adapting that result to our convention, the 
anomalous dimension of the gauge parameter is 

7a(/) = aoff^ + ai/ + 0(/) , 
1 /13 



= 16^ I^HT " "J " 3^^"^^ 
^ (lif (^^ (t " T" " ^"') " 10^^^/^^ - 8^^^^/^^) • (24) 
The anomalous dimension, 7^2(3^), of the composite operator Aj^ was calculated in |lj and reads 

7a2(<?^) = 709^ + 71/ + 0(9®) , 

70 = ^r^[{35 + 3a)CA-l6TFNf] , 

71 = ^ ^ [(449 + 33a + ISa^) Cl - ISOCATpNf - WlCpTpNA . (25) 

24 (167r2) LV / J 

In order to determine 6{g'^,a), we still require the counterterm 6^- In principle, this can be 
directly calculated from the divergences in W(J) when the propagator for a gluon with mass 
J is used. However, a less cumbersome way to compute 6( was described in It is based 
on the fact that the divergences arise in the 0{J'^) term and therefore that part of the Green's 
function which contains these divergences is equivalent to the Green's function with a double 
insertion of the operator. More specifically, one has two external J insertions with a 

non-zero momentum flowing into one insertion where the only internal couplings are those of 
the usual QCD action. Moreover, one does not require massive propagators but instead can 
use massless fields which simplifies the calculation. Therefore one is reduced to computing a 
massless two-point function for which the MiNCER algorithm, was designed. We used the 
version written in FORM, |47| I48j . where the Feynman diagrams are generated by Qgraf, |49j . 
to determine the divergence structure to three loops. Although we only require the result to two 
loops the extra loop evaluation in fact acts as a non-trivial check on the two loop result. This is 
because the emergence of the correct double and triple poles in e at three loops, in a way which 
is consistent with the renormalization group, verifies that the single and double poles of the two 
loop expression for 5C, are correct. We found 



+ -icA{-^^-la-\a'-laA+lTFN^]+0{g'), (26) 
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where A^^ is the dimension of the adjoint representation of the colour group. Assemblying our 
results leads to 



6{g ,a) 
So 



60 + 5ig^ + 0(/ 
Na 



167r2 
1 Na 



a 



6 (16^2y 



(Ca (-278 - 15q - 12q2 - Sa^) + 64rp7V/) . 



(27) 



As a check we see that d^g^^a) contains no poles for e ^ 0. Further, the expressions (|27() lead 
to the same results which were obtained earlier in the case of the Landau gauge (a = 0), as can 
be inferred from f9 without quarks and [T^ with quarks. 

From the renormalization group functions (ESJ, (HU), 1^5]) and (P7|). it is easy to see that the 
equation (|T6]) can be solved for by expanding C,{g'^,a) in a Laurent series as 



C(5',«) 



Co (a) 



+ Ci{a)+0{g^ 



(28) 



Substituting this expression in equation (fT?T|) . we obtain 

270 Co + So , 
270C1 + 271 Co + 5i • 
Thus (Eni) gives 



0/9 /- I '^Co 
2poCo + aao^ 



2/^*1 Co + aao-^ h aai^ 



dCo 



(29) 
(30) 



^ , . 2aCo + 3 (78 - 26a2 + Za^ + l8a In |a|) CaNa + 48 (a^ _ 3) NANfTp . 

Co(ct) = o 1 (31) 

^ ' 2((3Q-13)CA + 8iV/Tp)2 ^ ^ 

with Co a constant of integration. As a consequence of the already rather complicated structure 
of Co, we will determine Ci without quarks present corresponding to = since the expression 
for Ci with 7^ is several pages long. Using Mathematica, we find 



Ci(«) 



1 



12207367r2(13 - 3a)4 
+ 12716Coa2 I -442 - 132a + 540^ 



(-12207367r2a^^/^^ |-13 + 3a 

3a \ 4/13 



4/13^ 



1287(1 -- 



02^1 



4 4 17_3a 
13' 13' 13' 13 



Ca 



+ 



( 



1697175909 1 



3a \ 4/13 



13 / 



4 4 4 17 17_ 3a 
13' 13' 13' 13' 13' 13 



3335904 1 



3a \ 4/13 
13; 



17 17 17 30 30 3a 
13' 13' 13' 13' 13' 13 



+ 17 (-396870474 + 368850105a - 48761440a2 + 20662140^ + 1928718a4 
- 1004751a^ + 60588a*' - 12894024 



3a \ 4/13 2 ^ 



1-^) 

13 J 



9 9 4 3a 
i3'~13' 13' 13 



833976 



3a\4/i3 4 ^ r 4 17 30 3a 

1 a42Fi 

13/ 



13' 13' 13' 13 



3a \ 4/13 



3a\ 
13/ 



+ 2059992a^ In \a\ + 833976 
- 43758a^(-1961 + 7021n|a|)))iVA) , 



a42Fi 



8926632a^ ln|al 

^ 17 30 3a 
13' 13' 13' 13 



In jal 



(32) 
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where Ci is a constant of integration and the (generahzed) hypergeometric function is 

(ai)fc • • • («p)fc 



pFq[ai,--- ,ap;bi,-- - , hq] z\ = ^ 



(33) 



where 

{a)k = a{a + l)---{a + k-l) , (34) 

is the Pochhammer symbol. We note that Co (a = 0) = j^j^ and Ci(« = 0) = ^21^^ which 
recovers the Landau gauge results of [Sl^J. Further, the constants of integration Co and Ci do 
no enter the Landau gauge results. 



From expression ((T^. we deduce that the tree level gluon mass is provided by 

2 _ 9<y 



m 

while the 1-loop effective potential becomes 
2Co 



Co ' 



(35) 



1 - ) + -Indet 



Co^ 



1 



1 Cl 2 1 , 



2 

Na \ 
2 



^'^^(a^+m^) - 1 



a 



{d — l)tr In ( 9^ + ) + tr In ( 5^ + am 



(36) 



In dimensional regularization and using the MS scheme, one finds 



C: 



+ 



2C( 

64^2 ^2 



3Na g'a^ 



Co^ ; 647r2 Co' V 6 Co7^ 



-^ + ln.^^ 



■2 / ' 



Cor J 



(37) 



where ^ is the renormalization scale. It can be easily checked that the infinities in the effective 
potential cancel when the counterterms are included. 

Next, we look for a non-trivial minimum of the effective potential, which amounts to solving the 
gap equation ^ = 0. To avoid possibly large logarithms, we will set ^ = ""t-^ ~ Co ™ 
equation, 



dV 

da 



a 



"Co 



Co 



Ci 2\ , ^NAg^a / 5\ 3Na g^a 



Co- 



327r2 Co V ey 647r2 Co 



77 + 



+ 



Na a^g^a f 3 



32772 ^2 



-- + lna I + 



Na a'gV 

64772 ^2 



0, 



(38) 



and use the RGE to sum leading logarithms. Defining y = ^q^, we find as a solution of (|38|) 

CaCo 



o" = or y 



I6772C1 + ^ (1 + q2 - a2 ln|a| 



The first solution corresponds to the trivial vacuum, while the second one leads to 

3 

where the 1-loop formula for the coupling constant 



1 



/3oln^ 

MS 



(39) 



(40) 



(41) 
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was used. The vacuum energy is given by 



E = - 



1 Na 

2 647r2 



(3 + a^) 



4 

m . 



(42) 



We now consider the numerical evaluation of our results and restrict ourselves to the colour 
group SU{3). For SU{N) one has Tp = ^, Cp = Ca = N and Na = N'^ - 1. For 

completeness, we quote the results for the Landau gauge a = 0. 



y 



Landau 



.^Landau 
^Landau 



36 
187 

17 



0.193 



m = e24A5^Ri 2.031 Ams , 

^ re^Alj^!^ -0.323A: 



4 

MS 



(43) 
(44) 
(45) 



IQ^2- - ^ MS 

The results for general a are displayed in the Figures 1-3. For the moment, we have set 



0.15 
0.1 
0.05 




-1.5-1.25 -1 -0.75-0.5-0.25 



0.194 
0.192 
0.19 
0.188 
0.186 
0.184 




0.2 -0.1 0.1 0.2 



0.2 0.4 0.6 0. 



Figure 1: y as a function of a. 





-0.7-0.6-0.5-0.4-0.3-0.2-0.1 



-0.2 -0.1 



0.1 0.2 




Figure 2: m as a function of a. 




a 



Figure 3: ^^^ac as a function of a. 

Cq = Ci = 0. Evidently, y should certainly be positive and also relatively small to have a 
sensible expansion. Hence, we conclude from Figure 1 that we should restrict the range of 
values for a further. We also see that m becomes rapidly larger and E^^^ becomes rapidly 
more and more negative as a gets more negative. A more urgent problem is the fact that the 
vacuum energy E^.^^ depends on the gauge parameter a. Since E'vac is a physical quantity, it 
should be independent on the gauge parameter a. In the next section, we shall give a detailed 
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account of this gauge parameter dependence. We shall see that it is related to the impossibility 
of evaluating the effective potential to arbitrary high loop orders. Further, we shall provide a 
simple way to circumvent this problem and obtain a vacuum energy which is independent of a. 



3 Investigation of the gauge parameter dependence. 

One possible explanation as to why E^^^ depends on a could reside in the values of the constants 
of integration Co and C\ we have chosen. With another choice for these constants, it could be 
that E^^^ does not depend in a, or equivalently E^^^ = E^^^'^'^'^. This can be investigated by 
considering the expression H42() for E^^^{a, Co, Ci). In order to have the same E'^ac for each value 
of a, we should solve the following equation 



da da 







<^ a 



3 f 2^(dy 9y d(o dy d(i 
lly2 ' [da dCo da 8(1 da 







(46) 



in terms of Co and Ci. However, the solutions of this equation depend on a, and this is not 
allowed since Co and Ci should be a independent constants. This means that the a-dependence 
of i?vac cannot be eliminated by a suitable choice of Co and Ci . 

3.1 BRST symmetry and gauge parameter independence. 



Let us now turn to a more general analysis. Consider again the generating functional (|2()|) . We 
have the following identification, ignoring the overall normalization factors 



exp -iW{J) 



Dip] exp iS(j{J) 



1 



DipDcr] exp i 



S{J) + J d^x 



2C U 



42 




where S{J) and Sa{J) are given respectively by ©, and (|^T|) . Since 



d 1 
daN 



[Da] exp 



we find 



dWjJ) 
da 



d^x 



d^x 




-1 = 



cb 



da 



+ terms proportional to J 



which follows by noticing that 

dSjJ) 
da 



,4 fb^'b'' acj2' 

" X h t: 

I 2 da 2 , 



cb 



s d X [ — ) + terms proportional to J 



(47) 



(48) 



(49) 



(50) 



We see that the first term in the right hand side of (|50|1 is an exact BRST variation. As such, its 
vacuum expectation value vanishes. This is the usual argument to prove the gauge parameter 
independence in the BRST framework [2]. Of course, this is based on the assumption that the 
BRST symmetry is not broken. Notice therefore that there does not exist an operator Q with 
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= sQ, so that a non-vanishing vacuum expectation value for the condensate does not 

break the BRST invariance. Indeed, from 

9 
2 



sa = ^-sAl = -gAld^c- , (51) 



one can easily check that 

s j (fxC{A^,a) = , (52) 
so that we have a BRST invariant cj-action. 

The rest of the argument is based on the fact that J = when the vacuum is considered, so 
that we are left with only the BRST exact term in (|49() . More formally, the effective action 
V{a) = r is related to W(J) through a Legendre transformation 

rg)=->v(J)-/AJ(y)^. (53) 

The effective potential V{a) is then defined as 



a 



V{a) j d'x = T[-] . (54) 



Let amin be the solution of 



Since 



dV{(j) 
da 

5 



(55) 



r = - J , (56) 



one finds 



o- = amin ^ J = , (57) 
and invoking (|57j) . from (|53j) and (|54j) we derive 



da 



I d^x = ^W{J) 



j=o 



(58) 



Finally, combining (|49|) and (|58j) 

= . (59) 



From this, we conclude that the vacuum energy i^^ac should be independent of the gauge pa- 
rameter a. 

Apparently, our explicit result (|42|) for i^^ac is not in agreement with the above proof that i^^ac 
is the same for each a. If we examine the proof in more detail we notice that a key argument is 
that J becomes zero at the end of the calculation. In practice, this is achieved by solving the 
gap equation. Now, in a power series expansion in the coupling constant, the derivative of the 
effective potential with respect to a is something of the form 

(^^0 + vig"" + 0(/)) a , (60) 

where we assume that we work up to order g^ and that we have chosen JI so that the logarithms 
vanish. Then, the gap equation corresponding to (jHOJ reads 

vo + vig"" + 0(/) = . (61) 
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Due to (fKl|) and one also has 

J = g[vo + vig" + 0(/)) a . (62) 
This means that, if we solve the gap equation (|6T|) up to certain order, we have 

J = g{0 + 0(/)) a . (63) 

We also have 

da da g"^ da 
So, working to the order we are considering 

^ J = [^-0 + [-g^^vov, + — .0 j g +0{g))a . (65) 
From the square of the gap equation 

+ 2i;i^;o5' + 0(/) = , (66) 

it follows that 

j2 fdCl2 2 , rM.A\\ _2 



a^^'=lS"»''^°*''r^ (67) 

We see that, if one consistently works to the order we are considering, terms such as ^J"^ do 
not equal zero although J = to that order. Terms like those on the right hand side of H67() are 
cancelled by terms which are formally of higher order. This has its consequences for the terms 
proportional to J in ()49() . If one were able to work to infinite order, the problem would not 
arise. However, we do not have this ability, and we are faced with a gauge parameter dependence 
slipping into E^^^. 

3.2 Circumventing the gauge parameter dependence. 

We could resolve this issue by saying that the gauge parameter dependence of the vacuum energy 
should become less and less severe as we go to higher orders, and that eventually it will drop out 
if we go to infinite order. However, this is not very satisfactory, especially since we can surely 
never calculate the potential up to infinite order. Also as is clear from the quite complicated 
expression for Ci(q^)) which will enter the differential equation for C2(q^)> a 2-loop evaluation of 
the effective potential is already out of the question. 

Therefore, we could try to modify the LCO formalism in order to circumvent the gauge parameter 
dependence of -Evac- Therefore, we consider the following action 



S{J) = Sym + Sgf+fp + J df^x 



instead of (jEJ where, for the moment, JF(q) is an arbitrary function of a of the form 

J^{a) = 1 + Ua)g^ + h{a)g^ + Oig"") , (69) 
and J is now the source. The generating functional becomes 

exp -iW{J) = J [D(j)] exp iS{J) . (70) 
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Taking the functional derivative of VV( J) with respect to J, we obtain 



6W{J) 



SJ 



j=o 



-^(«)(- 



(71) 



Again, we insert unity via 



1 = - /[Z?cF]exp 



a 



^2 



to arrive at the following renormalized Lagrangian 



1 



From the generating functional 



1 



+ ^GF+FP 2g^J^{a)Z^C ^ 2 g^T{a)Z^C 



^2 _~ a2 1 ^2 



exp— iVV(J) = / [Z^^Jexpi 



it follows that 



5W{J) 



6J 



where the anomalous dimension of a equals 



7?(5) = ~^ = ^+7A^(5) + /^ 

The lowest order gluon mass is now provided by 

2 _ 90- 



9 In jF(a) 
3/i 



m 



Co ' 



and the vacuum configurations are now determined by solving 

dV{a) 



da 



. 



(72) 

-(73) 

(74) 
(75) 

(76) 

(77) 

(78) 



with V{a) the effective potential. In the MS scheme, the 1-loop effective potential reads 



2Co 



Ci 



Vi{a) = _ I 1- (2/o + ^ + - 



Co 



2 Na g"^ fS a 



£ 167r2 Co V 2 2 



3{NA)g^a^ 



+ 



2 5 gJ 
h In 



647r2 Co V 6 Co/^V 647r2 Co V 2 



+ 



2 3 ac/o- 
+ In —2 



+ 



a 

Na a'g^a ' 
647r2 Co 



1 - ( 2/o + ^ ) 5^ ) + 



mA9^ 

64^2 Co^ V 6 ' Co/I- 



Co 



3 , , aga 



^+ln.^^ 



(79) 



We included the counterterm contribution here to illustrate explicitly that Vi(ct) is finite. With 
((7B|) . it can also be checked that 



ti^V^{a) = ^ + 0{g^). 



(80) 
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Now, we can continue with the determination of the 1-loop vacuum energy, which will not only 
depend on a, Co and Ci, but also on /o(a). We will determine an expression for /o(a) so 
that E^^^{a,Co,Ci, fo{a)) does not depend on a. In the meantime, we could also absorb the 
constants of integration Cq and Ci in /o(a) so that E^^^ does not depend on them either. Based 
on this, we will immediately set Cq = Ci = 0. As usual, we put /I^ = ^ in the gap equation, 
which now reads 

dV 
da 



Co V V Cor y 327r2 Co^ V QJ 647r2 Co^ 

and use the RGE to sum the leading logarithms. One finds, in addition to the trivial solution 
a = 0, 

^ " 167r2(2/oCo + Ci) + ^(l + a2-a2in|a|) ' 

m = AmsB^ , (83) 

In principle, the analytic solution for /o(a) can be obtained by solving the following differential 
equation 

^ = ^ 2«m4 + 4(a2 + 3)m3^ = 
da da 

dy , dy dCo . dy dCi , dy dfo 



a 



/o , 2W ^y , <^y ^^0 , (Jy (71,1 , oy ujo \ , . 



lly2 c^Co da 9Ci c^/o da 



The quantity /o(a) constructed in this fashion will ensure i?vac(a) is independent of the gauge 
parameter a. However, we still have the freedom of choosing an initial condition. We will 
determine /o(a) so that E^^^^a) = i?vac(0) = -E^aT'^'"'- This amounts to choosing /o(a = 0) = 0. 
We can justify this choice based on our remark in the introduction, which is that Aj^ coincides 



with the gauge invariant quantity ^4^;^ in the Landau gauge in the FMR. Since our calculation 
is based on a perturbative expansion around A'^ = 0, which lies within the FMR, we stay within 

the FMR [sniisniiHi]- 

Unfortunately, the differential equation (|85|) is very hard to solve analytically. We could solve 
H85|) and consequently y, m and -Evac numerically. However, there is a more elegant way to 
obtain the analytical solution for fo{a). Considering the colour group SU{3) for simplicity, then 
since we know that by construction that E^^^ = E^^^^^'^'"'^ , we are able to write down the analytical 
solution for m as 

[s + a^) 

where use was made of and (|HH) . Putting (|5H)) in (|5S)) . we deduce that 



1/4 

m = ( — 5 1 Ams > (86) 



y = ^ • (87) 

661n3^ + 187 

Combining H82|) and ()87() finally gives the analytic expression for /o(a) 

^ (66 In + 18?) - 4 (1 + q2 - a2 In |a|) - 16712^1 
/o(a) = — ^ ^TT^T • 
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We have displayed /o(a), and m{a) for the range of values — ^ < a < ^ in Figures 

4-6. As a check, we have also plotted, in Figure 7, E^^^{a, fo{a)) as given in (|S11) to verify 
that Svac(o, /o(o)) = -E'vaT'^'"'- We observe several features. Firstly, although /o(a) has some 
singularities in — , the quantities y, m and E^^^ are completely regular functions of 

a. Secondly, the expansion parameter y remains relatively small, which makes our numerical 
predictions at least qualitatively trustworthy. Thirdly, we also see that the value for the tree 
level mass does not change spectacularly in the considered region. In the Feynman gauge q = 1, 
we have m'^^y'"""" = 1.89Att^. 



50 
25 

0- 
-25- 
-50- 
-75- 
-100 



Figure 4: /o as a function of a with — ^ < a < ^ . 
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Figure 5: y as a function of a with — ^ < a < 



Before ending this section, there are several other points. We have determined J-{a) with the 
renormalization scale Jl chosen in such a way that the logarithms vanish. Other choices of 'p. are 
of course also valid. We did not explicitly write this /I dependence of !F{a) in 



Also, the procedure we have described here applies of course at higher order. For example, at 
2-loops, /i(a) will be required to remove the a dependence. If we were to work to infinite order 
in , we could transform the action S{J) H68|) exactly into the action S{J) Q by means of the 
transformation 

J=^,- (89) 
The corresponding transformation for the a and a fields reads 

a = T{a)a , (90) 
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a 

Figure 6: m as a function of a with — ^ < a < 

-0 . 31i I 1 

-0.315 - 

I -0.32 - 
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-4 -2 2 4 

a 

Figure 7: i?vac as a function of a with — ^ < a < ^. 



which will transform the effective potential Voo(cr) exactly into Voo{cr). As such, the constructed 
vacuum energy will be the same in both cases and independent of the choice of a. 



4 Gluon propagator in linear covariant gauges. 

In the gluon propagator in the Landau was investigated, and a fit of the lattice results 
gave evidence for a gluon mass. In the Landau gauge, the lattice also gives evidence for the 
existence of a non-zero (^^^^ condensate, based on the discrepancy in the 10 GeV region, between 
the behaviour of the observed lattice gluon propagator and strong coupling constant and the 
expected perturbative behaviour. The results could be matched together using an operator 
product expansion analysis with a non-zero (^^^^ condensate [HIIZIIH]- A combined lattice fit 

resulted in {^^^^ ~ (1.64GeV)^. This quantity was obtained at a scale of 10 GeV in the 

2\ 



MOM renormalization scheme. Later, it was argued that this (A~,) condensate could be 

\ '^/OPE 

explained with instantons [Sj. 

One will notice that we did not give the estimate for (^^) itself. From the identification H22|l 
and using the relation and the explicit result (^U), one finds 

) = -^'^'^^ " - (0-29GeV)2 (91) 
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The extra minus sign arises because we have rotated from Minkowskian to Euchdean space time 
to make possible a comparison with the lattice. We used = 0.233GeV, which was the value 
obtained in We should be careful not to misinterpret the relatively big difference between 
(Af^) and (1911). Although our result is non-perturbative in nature, it is still obtained in 

\ ^/ OPE 

perturbation theory and as such it only gives information from the high energy region (or short 
range), while the OPE approach of Boucaud et al can only describe the low energy (or long 
range) content of (^^^^ It was already argued in |5] that (|^^^ can receive long and short 
range contributions. The minus sign in front of our result has to do with the regularization and 
renormalization of the quantity i-^n)- We refer to 0j for more details. 



To our knowledge, there has been little attention on the lattice to the gluon propagator in a 
general linear covariant gauge. Giusti et al managed to put the linear covariant gauge on the 
lattice (Sni IMl EHj • The tree level gluon propagator of Euclidean Yang- Mills theory with a 
linear covariant gauge fixing is given by 

^M-(9) = ^(V-(l-«)^) • (92) 
This can be decomposed into the transverse and longitudinal parts as 

q \ q J q q 

where D^lq"^) is q^ times the one used in [SUl |^ EH ES] . In general, one determines D^{q) via 
the projector 

Pj^Al) = Qf.qu . (94) 

If there is a tree level gluon mass m present, as in (f7H)). the Euclidean gluon propagator in linear 
covariant gauges reads 

D,uiq) = - (1 - ' (95) 

+ V q + otm^ J 

with the value of m given in H86|). The longitudinal part of this propagator is 



1 



-,4 



D'-iq) = P^MDAq) = -^—^ - (1 - ^ ) ' (96) 

D^{q) is plotted in Figure 8, again using A^^g = 0.233GeV. Of course, we should not attach 
any firm meaning to this plot, since we are only considering the tree level propagator and do 
not include any renormalization effects. If we could calculate the form factors, we would also 
inevitably encounter the problem of a diverging perturbation theory in the infrared region. We 
cannot make any conclusion about the behaviour of the propagator in the IR from the above. 
Many other (non-perturbative) effects can influence the propagators form in the IR. Nevertheless, 
it might be worth noticing that the longitudinal part D^{q) is not proportional to the gauge 
parameter. A similar behaviour was found by Giusti et al, see e.g. Figure 4 of |S21- This is 
already different from the perturbative prediction of massless Yang-Mills theory with a linear 
covariant gauge fixing |28j . 



5 Conclusion. 



We have considered Yang-Mills theories in linear covariant gauges and constructed a renor- 

V 



malizable effective potential by means of the local composite operator formalism for The 



17 



4 

3 . 5 

3 

2 . 5 

2 

1 . 5 

1 

. 5 



/ 

/ 



0.5 1 1.5 2 2.5 3 3.5 

Figure 8: Diiq) as a function of g with < q < 3.5GeV for a = (fat), a = 1 (dashed) and 
a = 4 (wide dashed). 



formation of the gluon condensate of mass dimension two is favoured since it lowers the vacuum 
energy. As a result, the gluons acquire a dynamical mass m. We discussed the gauge parameter 
dependence of the resultant vacuum energy and observed that this is due to the fact that we 
do not work up to infinite order precision, but have to truncate the perturbative expansion 
at a finite order. We explained how this gauge parameter dependence can be avoided by a 
modification of our method. 

Although there is limited lattice data available for the general linear covariant gauges compared 
with the Landau gauge, it would be interesting to calculate the form factor of the longitudinal 
and transverse part of the gluon propagator to make a more detailed comparison possible with 
the lattice results of [213 E21 ISnj • It would also be useful to have direct evidence from the 
lattice community that the ^^^^ condensate exists and that the gluons become massive, in 
analogy with the Landau gauge. A further point worth investigating is the possible existence 
of ghost condensates in the linear covariant gauges, as is the case in the Landau gauge [511 155] . 
These condensates can modify the gluon propagator further. 
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